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ABSTRACT 



O 

\ I discuss the properties of pions in "partially quenched" theories, i.e. those in 

^ ■ which the valence and sea quark masses, my and m5, are different. I point 

Q"^ ■ out that for lattice fermions which retain some chiral symmetry on the lattice. 



e.g. staggered fermions, the leading order prediction of the chiral expansion 
is that the mass of the pion depends only on my, and is independent of m^. 
Qh! This surprising result is shown to receive corrections from loop effects which 

are of relative size 7715 In my, and which thus diverge when the valence quark 
mass vanishes. Using partially quenched chiral perturbation theory, I calculate 
the full one-loop correction to the mass and decay constant of pions composed 
of two non-degenerate quarks, and suggest various combinations for which the 
\ prediction is independent of the unknown coefficients of the analytic terms in the 

chiral Lagrangian. These results can also be tested with Wilson fermions if one 
uses a non-perturbative definition of the quark mass. 



^Email: sharpe@phys.washington.edu 



1 Introduction 



This note is inspired by recent work of the SESAM collaboration, in which they have studied 
the light meson spectrum with two degenerate flavors of dynamical Wilson quarks []I| . They 
have calculated the masses of particles composed not only of sea quarks, but also of valence 
quarks with masses which can differ from those of the sea quarks. In the latter case they 
are studying a "partially quenched" theory. They find that, for their range of quark masses, 
can be well represented by a linear function of the valence and sea quark masses. For 
example, the mass of the non-singlet pion composed of two degenerate valence quarks of 
mass my, calculated with quark mass ms, takes the form 



Myy{mv,ms] 



= cymv + cstris , 

Wilson 



with cy ~ cs- The case my = ms corresponds to unquenched pions in two-flavor QCD, for 
which the pion mass vanishes in the chiral limit in the expected way. What is surprising 
about Eq. (|l|) is that if one works at fixed ms, but extrapolates my 0, then Myy does 
not vanish. Instead one must go to a negative valence quark mass, my = —{cs/cv)ms, to 
make the mass of the valence pion vanish. 

As explained in Ref. this peculiarity is easily understood in terms of the properties of 
Wilson fermions. The point is that chiral symmetry is completely broken by the lattice 
regularization, and so the value of hopping parameter k at which the valence quark mass 
vanishes depends on the parameters of the theory. In particular it depends on the sea quark 
mass, and so one should define a variable critical hopping parameter, Kc{ms). The quark 
masses in Eq. ([T|) are, however, defined as m = (2k)~^ — (2Kc(m5 = 0))~^, where udms = 0) 
is the critical value for which the unquenched pion mass vanishes. The negative quark masses 
are then an artifact of using ndnT'S = 0) instead of K,c{ms) in the definition of my.0 

My main aim in this note is to discuss how the results would change were one to use fermions 
for which some remnant of the continuum chiral symmetry survives discretization. What I 
have in mind are staggered and "domain- wall" fermions In the former case an axial 
subgroup of the SU{4:) chiral symmetry remains on the lattice, while in the latter the full 
chiral symmetry is broken only by exponentially small corrections. The only property of 
both types of fermion that I need is that these symmetries become exact when the lattice 
quark mass vanishes. If I then assume that the chiral symmetry associated with the valence 
quark is broken dynamically, with the formation of a non-zero condensate (qvQv), it follows 
that there will be a Goldstone pion whose mass vanishes when my = 0. In other words, the 
assumption of dynamical breaking of the valence quark chiral symmetry implies 

Myy(my = 0,m5) = 0. (2) 



^Even though one expects Kc to depend on ms, the strength of the dependence of found in Rcf. ^ is 
surprising, and has interesting impHcations for the extraction of physical quark masses 

■^For the sake of brevity, I will refer only to staggered fermions in the following, although all such references 
apply equally to domain- wall fermions. 
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Numerical evidence suggests that such symmetry breaking does occur for all values of the 
sea quark mass (including the quenched case ms oo). 

I now add to this the assumption of linearity, namely that Eq. (|Tp holds also for staggered 
fermions. These two assumptions then imply that cs = 0, so that 



Myy{mv,ms] 



= cymv ■ (3) 

staggered 



In other words, Myv is independent of the sea quark mass, at least for small ms where 
linearity applies^ If correct, this would be a surprising result. For example, one could obtain 
the mass dependence of the physical light pion mass without the need to work with physical 
sea quarks. This seems implausible on physical grounds. For one thing, the cloud of light 
mesons surrounding the pion depends on 777,5. But perhaps such effects are of higher order 
than linear in the expansion in the quark masses — indeed, loop effects in chiral perturbation 
theory lead to corrections to proportional to m'^ Inm^. 

In fact, I will show that such loop effects are enhanced in the partially quenched theory. 
Although the leading term does take the form of Eq. (j^), the dominant correction for small 
my is proportional to myms In Tny. For fixed ms this correction gets arbitrarily large relative 
to the leading order term as my — > 0. Thus, as one approaches this limit, the valence pion 
mass obtains a significant dependence on ms- This breakdown of Eq. occurs because the 
assumption of linearity fails, due to the appearance of non-analytic terms. The assumption 
that the valence pion mass vanishes when my —>■ remains valid. 

My main conclusion is thus that one cannot use Myy with unphysical sea quark masses to 
give an accurate estimate of the mass of the physical pion. What one can use Myy for, 
however, is to provide a sensitive test of loop effects predicted for the partially quenched 
theory. To this end, I have calculated the one-loop corrections to both the pion mass and 
decay constant, as a function of valence and sea quark masses, using partially quenched 
chiral perturbation theory |Q. In general these predictions depend on unknown constants 
multiplying the analytic terms of 0{m'^), but for certain combinations the analytic terms 
cancel. 

A similar deviation from linearity is predicted for Myy in fully quenched QCD. Using 
quenched chiral perturbation theory p, ^, one finds a correction proportional to r77y7r7Q In my. 
Here ttt-q is a parameter which, in the QCD chiral Lagrangian, gives the r]' its mass. There is 
some numerical evidence supporting this prediction, but the situation is muddied by the pos- 
sibility that finite volume errors mimic the chiral logarithm. For a review, see Ref. 0. The 
new predictions presented here provide another way of searching for chiral logarithms, and 
thus may help clarify the situation in quenched QCD. The advantage of partially quenched 
theories is that, for reasons explained in the following, the unknown parameter m^ does not 
appear in the predictions. 

This note is organized as follows. In Sec. ^ I give a brief description of the method of 
calculation, and then present the results for pion masses and decay constants in Sec. |^. 

^This result is for fixed bare coupling, go, because the "constant" cy depends upon qq. 



2 



Section ^ contains some general comments on predictions for partially quenched baryons, 
and Sec. ^ some conclusions. 

2 Calculation 

Partially quenched chiral perturbation theory has been described in detail in Ref . . I give 
here a summary of the aspects relevant to the present calculation. 

Consider a theory with two valence quarks, of mass mi and m2, and > 1 unquenched 
quarks of mass rris. To cancel internal loops containing the valence quarks one needs two 
ghost quarks (commuting quark fields q) with masses mi and 1712. For N = 2 this is the 
theory studied in Ref. Collecting all fields into a vector, 

Q = (qvi, qv2, qsi, qs2, qsn, gvi, qv2) (4) 

one sees that the chiral symmetry is the graded group SU{2 + N\2)l x SU{2 + N\2)r. The 
chiral Lagrangian consistent with this symmetry is 

C = ^str(a^Ea'^Et) + :^str(xSUSx)+"*'9^$o5^$o-m^<l>2, 

+ {«4 str str (xS^ + Sx) + 2^jia, str + Sx]) 

+a6 str (xSt + Sx) ' + ag str (xS^x^^ + Sx^x) } + .... (5) 

Here S = exp(22$//) contains all the Goldstone bosons, including the flavor singlet field^ 
$0 = str$/-\/3. The quark masses enter through x = 2/iM, where M is the mass matrix, 

M = diag(mi, m2, m^, m^, . . . , m^, mi, m2) , (6) 

with entries of ms. I have only kept those terms in the Lagrangian that will be required 
for the following calculations. In particular, I have omitted the arbitrary function of $0 
which can multiply each term. 

The terms multiplied by the coefficients a, are non-leading in the chiral expansion, since 
they contain an additional power of or m compared to the leading order terms.0 They 
give rise to the corrections to physical quantities which are analytic in the external momenta 
and quark masses. The other source of corrections is loop diagrams involving vertices and 
propagators coming from the leading order Lagrangian. These give rise to the non-analytic 
"chiral logarithms" , as well as analytic contributions. The typical size of both higher order 
corrections is p^/A^, where the chiral scale is A^ = Anf. 

^The factor of in <i>o is chosen so that, if the same normahzation were used in the corresponding QCD 
chiral Lagrangian, one would find ni^y — mQ/{l + «$) + 0{mq). 
^For a full list of these higher order terms see, for example, Ref. 
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The parameters in the chiral Lagrangian — /, /i, mo, a$ and the ctj — are not known a priori. 
They are functions of A^, the number of sea quarks, and thus are different for QCD {N = 3) 
and the N = 2 case considered in Ref. []I|. The most useful predictions of partially quenched 
chiral perturbation theory are thus those that depend on as few of these parameters as 
possible. 

The calculation of the one-loop corrections to the masses and decay constants is straightfor- 
ward. It is very similar at all stages to the corresponding calculation in quenched QCD, which 
has been discussed in Refs. ^. The only significant difference occurs in the propagators 
of mesons which have a flavor-singlet component. Consider, for example, the propagator of 
the meson with flavor composition qviQvi- In quenched QCD this is 

where Mf^ = 2/xmi is the leading order mass. The first term is the usual non-singlet prop- 
agator, while the second arises from "hairpin" diagrams in which the quark and antiquark 
annihilate. It is the second term which leads to enhanced chiral logarithms. The corre- 
sponding propagator in the partially quenched theory is different because of the possibility 
of internal loops of sea quarks. The result is [ffl 



p2 + Mf, 



(p2 + M2i)2 1 +jiV/3)K + «<i.p2)/(p' + Mis) 



p^ + Ml^ p^ + M\ 
1 



+ ~5 , ..2^J • (9) 



Here 



= ^ ' > ^ — ^ 10 
l + a$(A^/3) ^ ' 

is the mass of the singlet "77"' meson in the unquenched SU (N) sector of the theory. Similar 
results hold for the other flavor- singlet propagators. 

The r]' mass M does not vanish in the chiral limit. In the following, I will simplify the calcu- 
lation by assuming that the rj' is a heavy particle, i.e. M ^ A^, so that it can be integrated 
out of the theory. This is equivalent to assuming that the ratios Mg^/M^, Mf-^/M'^, etc. are 
small. This is certainly reasonable for N = 3, for which we know that M ^ M^/_phys ~ 1 GeV. 
Even for iV = 2, it is a sensible approximation, since the rj' will be comparable in mass to 
the vector mesons, which we do not include in the chiral Lagrangian. 

This assumption leads to two simplifications. First, integrals involving the rj' propagator 
+ M"^) can be dropped. These integrals can be expanded in powers of Mgg/M'^, 
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Mf^/M'^, etc. and their contributions can be absorbed by changing the coefficients in the 
chiral Lagrangian. For the quantities considered here the effect of r]' loops is to shift yU, ctg 
and ctg, as I have checked by exphcit calculation. Since we do not know these parameters 
a priori, we lose nothing by dropping the contribution from the f]' propagator. Indeed, this 
allows the results for iV = 3 to be matched directly onto those from the usual QCD chiral 
Lagrangian, from which the rj' has been integrated out. 

The second simplification is of the part of G^'^*^^ which remains after the rj' contribution has 
been removed. In this remainder, we can discard terms suppressed by powers of M^g/M"^, 
etc. for these are of the same size as two-loop terms which are not included. The propagator 
then simplifies to 

Note that the double pole term remains (and is the source of the enhanced chiral logarithms), 
but that the unknown parameters and a$ do not appear. In the unquenched theory 
{Mss = Mil), the propagator goes over to the usual form, with only a single pole, and 
with the term projecting against the r]'. The corresponding form for the off-diagonal 
propagator between a meson of composition qviQvi and qv2<lv2 is 



«S«V)«4f^^^ + |l"2,). ,12) 



3 Results 



I have calculated the complete one-loop correction to the mass and decay constant of a pion 
with composition qvi(lV2- I call these M12 and /12, respectively. Note that this state is a 
flavor non-singlet, and so there are no disconnected contributions to its propagator. Various 
limits of the general result are of interest: 



• If mi = m2 one obtains a non-singlet pion composed of degenerate valence quarks. I 
refer to the results in this limit as Mn and fu, or generically as Myv and fw 

• Setting m2 = rris one obtains a pion with only one quenched quark. I refer to the 
results for this pion as Mis and fis, or generically as Mys and fvs- 

• The case mi = m2 = mg is special, for then one is considering the physical unquenched 
pion in which the valence and sea quarks have the same mass.[] This case deserves a 
separate notation, and I follow Ref. [||] by denoting its mass and decay constant Mss 
and fssi respectively. 

^This is only true for iV > 2. One cannot make a non-singlet pion if there is only a single flavor of sea 
quark. Thus if = 1, the results for M12 in the limit that mi = m2 = ms do not correspond to those for 
an unquenched pion. 



5 



There is no problem in principle extending the calculation to pseudoscalars which have a 
flavor singlet component. I have not done so for two reasons. First, the results involve 
a number of unknown constants not present in the expressions for the non-singlet masses 
(an example is given in Ref. [^]). Second, nearly all simulations calculate only non-singlet 
masses, because the annihilation diagrams needed for singlet states require much greater 
computational resources. 

The leading order results for non- singlet pions are 
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tree 



/i(mi + 1712] 



12 



(13) 



The result for M12 agrees with that given in the Introduction, Eq. (^. In particular, M12 is, 
at leading order, independent of the sea quark mass m^. 



I have calculated the one-loop results using dimensional regularization and MS subtraction. 
To present these I use the notation that i/u = /i(mi + m2)/A^, yss = S/im^/A^, etc., where 
= 4:71 f. The result for the non-leading contribution to the pion mass is 



1— loop 



/i(mi + 7712) i ^ 



yujyss - yii) inyn - y22{yss - ^22) iny22 

2/22 - yn 



+ 2/12 (2^8 - as) + yssN{2a6 -a^)} . 



(14) 



The constants are to be evaluated at the scale Ay. In the degenerate limit this becomes 



= 2/imy I — [{2yvv - yss) In yvv + {yvv - yss)] 



+ yvvi^as - as) + yssN{2ae - a^)^ . 



(15) 



If we further set my = we obtain 



Mis 



i-ioop ^ ^'""^•^ |^2/ssln2/5S + yss [(Sag - as) + A^(2a6 - a^)] \ . (16) 



This agrees with the result from standard chiral perturbation theory ||^ . 

The enhanced chiral logarithmic corrections discussed in the Introduction appear in the result 
for Myy, Eq. (p!5|): these are corrections of relative size In my. Enhanced corrections are 
present also in the general result Eq. ([T^ ) if we take mi and m2 to zero in fixed ratio. They 
are absent, however, from the mass of the pion composed of one valence quark and one sea 
quark. This is obtained by setting mi = my and m2 = m^, yielding 



fi{my + ms) \ -T-jyvv In Vvv + yvsi'^ot^ - as) + yssN{2a^ - a^) 



l-loop ' ' ' ■ 

The chiral correction here is of relative size mylnmy, and thus vanishes when my —>■ 



(17) 
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The fact that all the loop terms are proportional to appears to be an accident. One 
might have expected terms proportional to A^, since there are N mesons of the form qvQs 
which can appear in loops. It turns out, however, that such contributions cancel in the final 
result. What remains is the contribution from loops involving "hairpin" vertices. 




The corresponding results for decay constants are 

^ = 1 - ^(2/15 In 2/15 + 2/25 In 2/25) 

'2/112/22-2/552/12, 2/11 , \ 

In h 2/12 - yss 

2/22 - 2/11 2/22 J 

+ -abyu + -a^Nyss , (18) 
= ^ ~ ^yvs^^yvs + ^a52/yv + ^a^Nyss , (19) 

■^^^ ' ^-(2/V5ln2/V5 + 2/5sln2/55) - 7^ (2/55ln^^ + 2/55 -2/vy J 



+ ^(^byvs + ^a-iNyss ■ (20) 

Comparing these with the results for masses, we see that the enhanced chiral logarithms 
survive here in fys but not in fw, which is the opposite of the situation for the masses. For 
both quantities the enhanced logarithms are multiplied by The decay constants do, 

however, have contributions proportional to A^, but from logarithms which are not enhanced. 



To give a sense of the size and form of the corrections, I display the results for the VV, VS 
and SS pions for N = 2. To convert meson masses to physical units I take / = 0.1 GeV 
(which is the approximate value for this constant in QCD). To convert quark masses to units 
of the physical strange quark mass, rrist, I assume the leading order result /imst = ^l-,phys- 
This ignores the shift due to the one-loop corrections, but, as we will see, these corrections 
are of moderate size. Finally, I set the analytic constants, a^s, to zero. This is the simplest 
choice given that we do not know what values to use for N = 2. I have checked that the 
essential features of the plots are unchanged if are set to the values they take in QCD. 

Figure |1] shows the one-loop predictions for the masses of the VV and VS pions plotted 
against the valence quark mass my, for three values of the sea quark mass, = mst/4, 
mst/2 and mgt. This range is chosen to cover the typical values for "light" quarks used in 
present simulations. The chiral expansion is likely to break down at the upper end of this 
range. By construction, the VV and ^5* curves must cross when my = ms- For purposes 
of comparison, I have also included the one-loop result for plotted against ms- 

Lowest order chiral perturbation theory predicts that all curves are linear, with the three 
for Myy and that for M55 coinciding, while the curves for Mys have half the slope of those 
for Myy. We see that, although one-loop corrections do change this prediction, the major 
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Figure 1: Predictions for pion masses using values for the parameters discussed in the text. 
The sohd and short-dashed curves show Myy and My^, including one-loop contributions, 
plotted against my. The three sets of curves correspond to nis = nist, mst/2, and mst/4 
as one moves from top to bottom. The long-dashed curve is the result for M^g at one-loop 
plotted against m^. 

features of the leading order result remain. The most significant change is that Myv and 
Mss no longer coincide, with the curves for Myv showing some curvature. 

To magnify the difference between Myv and Mss, it is advantageous to consider quantities 
in which the leading order quark mass dependence has been removed. One such quantity, 
ln[M^2/ f^i^i + "^2)]; is plotted in Fig. The enhanced chiral logarithm causes the VV 
curves to diverge as my — > 0. A similar divergence is predicted for quenched QCD, and 
plots of this kind are a useful way of searching for this divergence. Such a search will 
not, however, be easy. For one thing, it is hard to distinguish logarithms from the linear 
dependence predicted by analytic terms unless one has a large range of quark masses.^ It is 
also difficult to separate logarithms from the 1/my behavior associated with finite volume 

*One cannot overcome this problem by considering the difference between the VV and SS curves — this 
cancels the analytic term proportional to 2as — as, but leaves the term proportional to 2aQ — a^. 
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Figure 2: Predictions for \nMyy/2fj,mv, In My^//i(my + rris) and \n Mgg/2^ms, including 
one-loop contributions, plotted against the average quark mass. Notation as in Fig. |T[ The 
three sets of curves for VV and VS mesons correspond to ms = mst/4, mst/2 and rrist as 
one moves from left to right. 



effects [g g. 

An alternative approach is to consider quantities in which the analytic terms cancel. One 
example is the difference between the masses of VV and VS mesons composed of quarks 
having the same average mass evaluated at the same sea quark mass. This vanishes at tree 
level, but not at one-loop: 

Mi\(mi, ms) - Mls{m2 = 2mi - rris, ms) = 2/imi^ 1 7/22 In— + 7/22 - yii \ ■ (21) 

t Z/22 J 

This may also be a more practical difference to study in detail as it relies only on changing 
the valence quark mass. Once one has determined /, this difference is, for small enough 
quark mass, a prediction of partially quenched chiral perturbation theory free of unknown 
parameters. To illustrate this prediction I have included the results for the VS mesons 
in Fig. ^ but now plotted against the average quark mass {my + ms)/2 rather than my- 
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The quantity in Eq. ( [2T| ) is simply the difference between the VV and VS curves at fixed 
rris- This brings out a striking (and seemingly accidental) prediction of chiral perturbation 
theory: the curves for Myy and My^ have the same derivative at my = ms- This is true at 
leading order and is not affected by loop corrections. Because of this, the predicted difference 
between the VV and VS curves is small (and gets smaller as the quark masses decrease). 

One can also consider a generalization of this difference which allows more flexibility in 
testing the predictions of chiral perturbation theory. Imagine working at fixed ms, and pick 
three valence quark masses which satisfy 

mi + m^ = 2m2 . (22) 

The quantity of interest is the relative difference between the masses of the "13" and "22" 
mesons. This vanishes at tree-level, but at one-loop takes the form 

- _ 1 j yujyss -yn)Hyn/y22) yu-yss , 

Mf3 + M|, -2iV\ y,,-yn 2 + ^ j ' ^^^^ 

In evaluating the r.h.s. of this expression it is legitimate to use the pion masses themselves, 
rather than their lowest order expression in terms of quark masses, since the difference is a 
two- loop effect. The same is true for Eq. (^). 

The results for decay constants are shown in Fig. ^. I have plotted them against the average 
quark mass so that the analytic contributions cancel in the difference between fw and 
fvs at fixed ms.[| It turns out that, as for pion masses, the VV and VS curves have the 
same derivative at my = ms- Note that the corrections are of moderate size even for 
mv,ms ~ mst. The only exception is that fvs diverges in the limit my due to the 
enhanced chiral logarithm. 

One way of looking for this enhanced logarithm is to form the ratio 

Rbg = (24) 

V ill 722 

at fixed ms- This quantity was introduced by Bernard and Golterman as a way of testing 
chiral perturbation theory for quenched QCD 0, but turns out also to be useful in partially 
quenched theories. The analytic contributions cancel in Rbg^ and one finds 



Rbg — 1 




1/111/22 , yn , 

In h yi2 



1/22 - yil 1/22 



- yss 



1/12 j^yn^^ 



1/22 - 1/11 1/22 



(25) 



It is again legitimate to use the actual pion masses when evaluating this result. Note that the 
prediction for this quantity diverges as mi with fixed m2 and m^. A similar divergence 
is predicted for quenched QCD, with ygg/2N replaced by mp/SA^ p. 

^It should be borne in mind, however, that the detailed shape of the curves, and the difference between 
the VV and SS results, does depend on the choice of 04 and as. 
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(mi+ni2)/2m^trange 

Figure 3: Results for decay constants plotted against average quark mass. Notation as in 
Fig. ^ As one moves from top to bottom, the three sets of curves correspond to = nist, 
mst/2, and mst/4. 



4 Baryon masses 



The enhanced chiral loops which lead to the divergences also appear in other quantities. In 
particular, one can use partially quenched chiral perturbation theory to study the behavior 
of baryon masses, using a straightforward extension of the methods developed for quenched 
baryons |TT|]. I have not carried out a detailed calculation, but it is easy to determine the 



general form of the dependence on the quark masses. The result for baryons composed of 
three degenerate valence quarks is 

Mvvv = Mo + CiMvvMls + c^M^y + csM^^ + O^M^y) . (26) 

This form applies for both spin 1/2 and 3/2 baryons, although the coefficients (including 
Mo) depend on the spin. The coefficient ci can in principle be predicted in terms of the 
pion-nucleon couplings F and Z), and the decay constant /. The other coefficients are 
A^— dependent unknown constants. The term proportional to Ci is the analogue of the 
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enhanced chiral logarithms found above. It does not diverge as my 0, but it is the 
dominant correction for sufficiently small my at fixed ms- This is no longer true in the 
unquenched limit, my = nis, for then the Ci term is of 0(M|^). 

5 Conclusions 

Partially quenched theories are a step on the way from quenched to full QCD. They allow one 
to partially probe the dynamics of light quarks by sending the valence quark mass towards 
zero, while holding the sea quark mass fixed. In this paper I have investigated the errors that 
this procedure introduces. The situation turns out be rather subtle in that the errors only 
show up at non- leading order, but they nevertheless diverge (in relative size) as my — *• 0. I 
have suggested a number of combinations of pion masses and decay constants with which to 
search for such divergences. 

The entire discussion has assumed that we know what the lattice quark masses are, as is 
the case for staggered fermions. As mentioned in the Introduction, with Wilson fermions 
there are problems in determining the quark mass from the hopping parameter in partially 
quenched theories. These problems can be avoided, however, by determining the quark mass 
non-perturbatively using the PCAC equation. In this way the predictions can be tested 
for Wilson fermions. The only disadvantage compared to staggered fermions is that the 
predictions will hold up to discretization errors of 0{a) rather than O(a^). Even without 
a non-perturbative determination of the quark masses, Eqs. (pT]), (p3D and ( pS]) can still be 
tested with Wilson fermions. 
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